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Abstract 

Effects of time-dependent applied fields on quantum coupled double-well (DW) systems with 
Razavy’s hyperbolic potential have been studied. By solving the Schrodinger equation for the 
DW system, we have obtained time-dependent occupation probabilities of the eigenstates, from 
which expectation values of positions aq and X 2 of particles ((aq + aq)), the correlation (r(t)) 
and the concurrence ( C(t )) expressing a degree of the entanglement of the coupled DW system, 
are obtained. Analytical expressions for (xi + aq), T(t) and C(t) are derived with the use of the 
rotating-wave approximation (RWA) for sinusoidal fields. Model calculations have indicated that 
(aq + aq), T(f) and C(t ) show very complicated time dependences. Results of the RWA are in good 
agreement with exact ones evaluated by numerical methods for cases of weak couplings and small 
applied fields in the near-resonant condition. Applications of our method to step fields are also 
studied. 
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I. INTRODUCTION 


Extensive studies have been made for quantum double-well (DW) systems in physics and 
chemistry where a tunneling is one of intrigue quantum phenomena [T]. Effects of applied 
fields on DW systems have been studied (for review see [2]). Various phenomena such as 
a coherent destruction of tunneling by applied fields were pointed out [3]. The two-level 
(TL) system which is a simplified model of a DW system, has been employed for a study 
on qubits which play important roles in quantum information and quantum computation. 
Many theoretical studies on effects of fields applied to single and coupled qubits have been 
reported with the use of the TL model HH3- In contrast to the simplified TL model, studies 
on coupled DW systems which are commonly described by the quartic potentials are scanty 
|8j, because a calculation of such a system is much tedious than that of the coupled TL 
model, even for the absence of applied fields. One of difficulties in studying coupled DW 
systems is that one cannot obtain exact eigenvalues and eigenfunctions of the Schrodinger 
equation for quartic DW potential. Then one has to apply various approximate approaches 
such as perturbation and spectral methods to quartic DW models. Razavy [2] proposed 
quasi-exactly solvable hyperbolic DW potential for which one may exactly determine a part 
of whole eigenvalues and eigenfunctions. A family of quasi-exactly solvable potentials has 
been investigated mm- 

Recently the present author [12] has investigated the relation between the entanglement 
and the speed of evolution in coupled DW system described by Razavy’s potential. It 
would be interesting to study effects of applied fields on coupled DW systems with Razavy’s 
potential, which is the purpose of the present study. Some sophisticated methods like the 
Floquet approach have been developed in solving Schrodinger equation for time-dependent 
periodic fields. In order to treat the periodic as well as non-periodic dynamical fields, we solve 
in this study the time-dependent Schrodinger equation by a straightforward method. An 
advantage of our approach is that we may exactly determine eigenvalues and eigenfunctions 
of driven coupled DW systems. We calculate expectation values of various quantities such 
as positions of particles, the correlation and the concurrence which expresses a measure of 
the entanglement of a coupled DW system. Effects of applied fields are analytically studied 
with the use of the rotating-wave approximation (RWA) which has been widely adopted 
for sinusoidal periodic held, in particular for the TL model. The validity of the RWA may 
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be examined by a comparison between results of the RWA and exact ones evaluated by 
numerical methods. 

The paper is organized as follows. In Sec. II, the calculation method employed in our 
study is explained with a brief review on Razavy’s hyperbolic potential [9j. Equations of mo¬ 
tion for populations of four energy levels are obtained from the time-dependent Schrodinger 
equation of driven coupled DW systems. Expressions for expectation values of particle po¬ 
sitions, the correlation and the concurrence are calculated. For sinusoidal fields, we present 
their analytical expressions by using the RWA. In Sec. Ill, we report model calculations 
with the use of the RWA and numerical methods when the sinusoidal fields are applied to 
the initial ground state. In Sec. IV, calculations are made for sinusoidal fields applied to 
the initially wavepacket state. Our method is applied also to the case of applied step fields. 
Sec. V is devoted to our conclusion. 


II. COUPLED DOUBLE-WELL SYSTEM WITH RAZAVY’S POTENTIAL 


Calculation method 


We consider coupled two DW systems whose Hamiltonian is given by 


// = //„ + //, 


c 


Hi, 


( 1 ) 


where 


h„ = y. 

n= 1 L 

H c = - gxix 2 


h 2 d 2 


2m dx\ 


+ V{x n ) 


Hi = -Oi +x 2 )F(t), 

e 


h 2 

v <*> - ^ 


e 


— cosh Ax — 4£ cosh 2x - 

o <> g 


( 2 ) 

( 3 ) 

( 4 ) 

( 5 ) 


Here X\ and x 2 stand for coordinates of two distinguishable particles of mass m, H 0 signifies 
a DW system with Razavy’s potential V(x) [9], He means the coupling term with an inter¬ 
action g , and Hi includes the time-dependent applied field F(t) whose explicit form will be 
given shortly [Eq. (40) or (88)]. The case of Hi which is more general than Eq. @ will be 
studied in the Appendix. The potential V(x) with h = m = £ = 1.0 adopted in this study 
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is plotted in Fig. 1(a). Minima of V{x) locate at x s = ±1.38433 with V(x s ) 
its maximum is V(x) = —2.0 at x — 0.0. 

Firstly we consider only H 0 in Eq. whose eigenvalues are given [S] 


eo 


ei 


£2 


^3 


2 m 
2 m 
2 m 
2 m 


-e-5-2v/4-2e±e 
e - s - 2 V / 4±2e±e], 
-e-5-2/f-2£-e 
e-5 + 2v'4 + 2e + e 2 l , 


5 


5 


and whose eigenfunctions are given by 


-8.125 and 


( 6 ) 

(7) 

( 8 ) 
(9) 


0 o(x) 

M x ) 

M x ) 

M x ) 


H 0 e -« cosh2a; / 4 cosh x + (4 - f ± 2 a/ 4 — 2£ ± £ 2 ) cosh 3x 

Ai e -« cosh2a; / 4 sinh x ± (4 ± f ± 2^4 + 2£ ± £ 2 ) sinli 3x 

A 2 e S c ° sh2x / 4 C osh x ± (4 - f - 2 a/4 — 2 £ ± £ 2 ) cosh 3x 

H 3 e ~e*»h2x/4 [ 3 ^ sinh x ± (4 ± f - 2^4 ± 2 £ ± £ 2 ) sinh 3 x 


( 10 ) 

( 11 ) 

( 12 ) 

(13) 


A n (n = 0,1) denoting normalization factors. Eigenvalues for the adopted parameters are 
eo = —4.73205, ei = —4.64575, 62 = —1.26795 and e 3 = 0.645751. Both eo and e 3 locate 
below V(0) as shown by dashed curves in Fig. 1(a), and 62 and e 3 are far above e±. In 
this study, we take into account only the lowest two states with eo and ei, which is justified 
because of e\ — e 0 (= 0.0863) <C e 2 — e 3 (= 3.3778). Figure 1(b) shows eigenfunctions of 
0 o ( x ) and 0 i (x), which are symmetric and anti-symmetric, respectively, with respect to the 
origin. 

Secondly we include the coupling term He in Eq- (|3]) . With basis states of 0o(xi)0o(x 2 ), 
0o(xi)0i(x2), 0i(aq)0o(^2) and 0i(xi)0i(x 2 ), the energy matrix of the Hamiltonian of H 0 ± 
He is expressed by 


Ho + Uc 


7 2eo 0 0 — gj 2 ^ 

0 e 0 ± ei —#7 2 0 

0 -g'f 2 e 0 + ei 0 

V ~ 91 2 0 0 2ei / 


(14) 


with 


7 = 


0o(x) x 0i (x) dx = 1.13823. 


(15) 
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FIG. 1: (Color online) (a) Razavy’s DW potential V(x) (solid curve), dashed and chain curves 
expressing eigenvalues of eo and ei, respectively, for H = m = £ = 1.0 [Eq.([5])]. (b) Eigenfunctions 
of 4>o{x) (solid curve) and <j> i(x) (dashed curve). 

Eigenvalues of the energy matrix of hio + hie are given by 

E 0 = e - V^ 2 + g 2 Yi (16) 

Ei = e- 37 2 , (17) 

E 2 = e + 37 2 , (18) 

E 3 = e + \/5 2 + 2 7 4 , (19) 

where 


e = a + e 0 = -9.3778, 
<5 = €l — e 0 = 0.0863. 


Corresponding eigenfunctions are given by 


$0(27, ^2) 

$i(xi,x 2 ) 

$ 2 (^ 1 , X 2 ) 

$ 3 ( 27 , x 2 ) 


cos 9 0 o (^i)0o(x 2 ) + sin# 0i(xi)0i(x 2 ), 
-j= [M x i)M x 2) + M x i)M x i)}, 

-J= [~M X i)M X 2) + M X l)M X 2)} , 

-sin# (j) 0 ( 27 ) 00 ( 27 ) + cos# 0i(xi)0i(x 2 ), 


where 


tan 2# 


9± , 

6 ' 1 


7T 

4 


<s< J) 


( 20 ) 

( 21 ) 


( 22 ) 

(23) 

(24) 

(25) 


(26) 
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We hereafter assume g > 0 [12]. The g dependence of E v (y = 0 — 3) is shown in Fig. 2 of 
Ref. P2j- 

Thirdly we take into account Hj for an applied held in Eq. Q. The energy matrix of the 
time-dependent total Hamiltonian H (= H 0 + He + Hj) with basis states of (fto(xi)(fto(x 2 ), 
4>o(xi)4>i(x 2 ), 4 >i(xi)(j)o(x 2 ) and is expressed by 


Hq + He T Hj 


( 2 e 0 - 7 F(t) —(? 7 2 ^ 

- 7 F (t) e 0 + ei -gj 2 - 7 F(t) 
- 7 F(t) -g~i 2 e 0 + ei - 7 F(t) 

V ~91 2 ~lF(t) ~lF(t) 2ei / 


(27) 


Alternatively, the energy matrix of H may be expressed with basis states of $ 0 , ‘hi, $2 


and $3 in Eqs. (22)-(25) by 


where 


H 


/ 

E 0 

— aF(t ) 

0 

0 


—aF(t ) 

E\ 

0 

~PF{t) 


0 

0 

e 2 

0 

V 

0 

~PF(t) 

0 

E 3 


a = v / 27 (cos 6 l + sind), 
fd = \/ 27 (cos 6 ) — sind). 


(28) 


(29) 

(30) 


In onr following analysis, we adopt the energy matrix given by Eq. (J28|) because it has more 
transparent physical meaning than Eq. 


We expand the eigenstate '&(x 1 ,x 2 ,t) of H in 
terms of $ v {x\,x 2 ) {y = 0 — 3) with the time-dependent expansion coefficients a v (t) as 

3 

T(f) = T(xi,x 2 ,t) = ^2 a v (t) § u (x ll x 2 ) (31) 


v=0 


where expansion coefficients satisfy the relation 


£M‘)I 2 = i. 


(32) 


v =0 


The Schrodinger equation: ihd^(t)/dt = H^(t) becomes 

3 , 3 3 

ih Y,J t = Y / a At)K^e- iEulh + J2 a At) H i^ e ^ iE, ' t/h - (33) 

u=0 L v=0 u=0 
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Multiplying <f>* (// = 0 — 3) from the left side of Eq. (33) and integrating it over x\ and X 2 , 
we obtain equations of motion for a M (t) 


mFFl 

dt 






(Xy (t) , 


(// = 0 - 3) 


(34) 


u =0 


where 


E v 


h 


(35) 


With the use of the energy matrix in Eq. (28), equations of motion for a M (t) become [the 
argument t in a u (t) is hereafter suppressed] 


ih = — aF(t ) e lAlot ai, 

ih^- = —aF(t) e iAlot ao — /3F(t)e~ iA3lt 
^ d a 2 

ih —— = o, 

dt 

ift—j ~ = —fiF{t)e lA3lt ai. 


0*3, 


(36) 

(37) 

(38) 

(39) 


When we apply the sinusoidal held given by 


F(t) = f sin uit, 


(40) 


Eqs. (|36|)-(|39|) become 

aQ [ e ^- Ai oh _ e -i( w +A 10 )t] ( 41 ) 

0 [ e -(-+ A io)i _ e -i(a,-A 10 )t] ao + [ e h-A 31 )t _ e —i(u+A 31 )t] a3> (42) 

(43) 


da q / ctj 

dt \ 2h y 

da i 
dt 


da,2 _ 

"rfT ~ ’ 


^ [ p i(uj+A 31 )t _ p — i(w—A 3 1 

7 L 


ai, 


dt \2h / 

where / and a; denote magnitude and frequency, respectively, of the applied held. 


(44) 


Rotating-wave approximation (RWA) 

In the rotating-wave approximation (RWA) where only terms with (a; — Aio) are taken 
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into account in Eqs. (41)-(44), we obtain 


dt - { 2h 1 e ai ’ 


dai 
dt 

da 2 c /03 




= 0. 


(45) 

(46) 

(47) 


dt dt 

For a given initial condition of a v (t) = a v ( 0) at t — 0 [y = 0 — 3), we obtain the solution of 


Eqs. (45)-(47) 


ao(t) = e *(^- A io)i / 2 _ Al0 _ e <nt /2 + _ Aiq + e -<nt/ 2 ] , ( 48 ) 


- (§) •*- 

Oi(f) = e -i(»-A.0)t/2 ( j4e fflt/2 + Be - i ni/2) j 


a 2 (l) = a 2 (0), 
a 3 (t) = a 3 (0), 


(49) 

(50) 

(51) 


with 


A = 

B = - 


m) ao(0) + ( 2 ^) (w ~ Ai ° +n) “ i(0) ’ 

^)“»(°)-( ^(-^-njaifO). 


where fl stands for Rabi’s frequency given by 


- A 10 ) 2 + ( af/hy. 


For a later purpose, we may rewrite a 0 (t ) and Gq(f) as 


a 0 (t) = rQ e ^- Al o-W + Sq e ^-Aio+a)t/ 2 , 


a 1 


— ri e ~i(u-A 10 -n)t/2 _|_ Sl e -i(u-A 10 +n)t/2^ 


(52) 

(53) 


(54) 


(55) 

(56) 


with 


r o — B ^ ( w — ^10 + 

^( 57 ) (u,_Aio “ !J) ’ 


s o 

ri 


= A, Si = B, 


(57) 

(58) 

(59) 


where r 0 , s 0 , rj and Si are time independent. 









B. Various physical quantities 


Once time-dependent a u (t) are obtained from Eqs. (36)-(39) or from Eqs. (55) and (56), 
we may evaluate various physical quantities such as expectation values, the correlation and 
concurrence. 

(1) Expectation values 

Time-dependent expectation values of (aq) and (x 2 ) are expressed by 


{x i) = 


'•OO /» OO 


' —OO J — OO 


+ (cos 6 — sin #)Re{aga 2 e 


—iA2ot 


, x 2 ,t) dx i dx 2 , 

(60) 

,—iAioi _ a * a3e -i^32ty 

+ a\a 3 e * A31 *}], 

(61) 

~ iAlot + a* 2 a 3 e~ iA32t } 

+ a\a 3 e * A3lt }], 

(62) 

e -iAiot| 

(63) 


Substituting Eqs. (55) and (56) to Eq. (63), the expectation value in the RWA with 


a 2 (0) = a 3 (0) = 0 is given by 

(xi )rwa = (x 2 )rwa = v^ylcos 9 + sin6>) Re{(r*si + sjjjri) e 
+ r*rq e-^-W + s* 0 s 1 e - i{uj+n)t }, 


—iut 


(64) 


which includes time-dependent components with frequencies of caiO besides u of the applied 
field. 

(2) Correlation 

The correlation T(t) is defined by [12] 

/ OO /»oo 

/ ^*(x 1 ,x 2 ,0) ^(x 1 ,x 2 ,t) clxi dx 2 \ 2 , (65) 

-OO J — OO 

3 

= | aS(0)a 0 (t) + J] al(0)a u (t) e~ lA ^ | 2 , (66) 

V =1 

which is unity at t — 0. 

In the RWA, the correlation is given by 

fmW 2 = l a o(0)| 2 |a 0 (t)|- + |ai(0)|-|ai(t)| 2 + 2Re{a 0 (0)aS(t)aJ(0)ai(t) e~ lAlot }. (67) 
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With the use of Eqs. (55) and (56), the correlation in the RWA is expressed by 


r RWAit) 2 — (?o + s o)( r o + s o)(l r o| 2 + |So| 2 + 2Re{sQroe ifM }) 

+ ( r i + s i)( r i + si)(|n| 2 + |si| 2 + 2Re{rj|'sie _int }) 

+ 2Re{(r 0 + So)(r^ + s^)[(sQr! + TqSi) e lujt 

+ r^ne^-^+ s* 0 s 1 e~ i ^ +Q)t ]}, ( 68 ) 


which consists of components with frequencies of oj, hi and cn±fl When we take the average 
of T RWA(t) 2 over a long period, oscillating terms vanish and its average becomes 


(r RWA(t) 2 )av = lilTlT^-oo / ^RWA(t ) 2 dt, 


(69) 


= ( r o + s o)( r o + s o)(l r o| 2 + l s o| 2 ) + ( r i + s i)( r i + -Si)(|ri| 2 + |Si| 2 ). (70) 


(3) Concurrence 


Substituting Eqs. (22)-(25) into Eq. (31), we obtain 


1^) — cqo|0 0) + cqi 10 1) + Cio11 0) + cn|l 1), 


(71) 


with 


coo 

Coi 

ClO 

Cll 


= a 0 cos 6 e lEot — a 3 sin 6 e lEat , 

= —j=( ai e~ iElt — a 2 e~ iE2t ), 

V2 


1 


— 7 =(ai e~ iElt + a 2 e ~ iE2f ) 


= a 0 sin 6 e lEot + a 3 cos 9e iEat , 


(72) 

(73) 

(74) 

(75) 


where |k £) = 4 > k{ x \)4 > (.( x 2) with k,£ — 0,1. The concurrence C of the state |T) given by 


Eq. (71) is defined by 


C 2 — 4 IcooCn — CqiCio| 2 . 


(76) 


The state given by Eq. becomes factorizable if and only if the relation: c 0 oCn—c 0 iCi 0 = 0 


holds. Substituting Eqs. (|72h-( 75) into Eq. (76), we obtain the concurrence 


C(t ) 2 = \[a 0 (t) 2 - a 3 (t) 2 e ~ 2iA30t ] sin 29 + 2a 0 (t)a 3 (t) cos 26 e~ iA:iot - a^t) 2 e~ 2iAwt 

+ a 2 (f) 2 e“ 2iA2(,t | 2 . (77) 
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In the RWA with 02 ( 0 ) = 03 ( 0 ) = 0, the concurrence is given by 


C'RWAit ) = | a 0 (t) sin 26 — ai(t) e 


2 -2iAi 0 O2 


( 78 ) 


With the use of Eqs. (55) and (56), Eq. (78) becomes 


Crwa ^) 2 = Dq — Di sin 2 d + D 2 sin 2 26, 


(79) 


with 


D 0 

D l 


N 4 + |si | 4 + 4|n | 2 | Sl | 2 + 2Re{2(|n | 2 + |si| 2 )rj(si e~ iQt + rf s\ e” 2 ^}, (80) 

2Re{ (r* 2 s 2 + s^r 2 + 4r*s*riSi) e~ 2iuit + 2 (r*s* 0 rl + rfr^) er^-W 


+ 2 (rQSoSl + e 


1 *Z Z 

+ r o r \ e 


I *z z 

+ Sg e 




D 2 = |r „| 4 + |s „| 4 + 4|r 0 nsor + 2 Re{ 2 (|r „| 4 + |»„|>;r„ e~ tm + f^r 4 e"““}, 


(81) 

(82) 


which include contributions from multiple components with frequencies of ui, 2 ui, fl, 2 h 2 , 
u ± 12, 2u ± 12 and 2 (cj ± 12). The concurrence averaged over a long period is given by 

(C RWA (t) 2 ) av = (|r 1 | 4 +| Sl | 4 + 4|r 1 | 2 | Sl | 2 )+sin 2 2d(|ro| 4 + | S o| 4 + 4|ro| 2 | So | 2 ). (83) 


III. MODEL CALCULATIONS 

Assuming the initial ground state given by 

a 0 (0) = 1, ai(0) = a 2 (0) = a 3 (0) = 0, (84) 

we have made numerical calculations by changing model parameters of f, g and c 0 . 


A. / dependence 


Figures | 2 ]( a), (b) and (c) show time developments of populations of |a„(t )| 2 in levels v 
(= 0,1,3) for / = 0.0, 0.01 and 0.02, respectively with g = 0.01 and oj = A 10 (=0.07431) 


obtained by numerically solving Eqs. (41)-(44) which is hereafter referred to as an exact 


calculation: note that 02 ( 2 ) = 0 because of d a 2 {t)/dt = 0 and 02 ( 0 ) = 0. For comparison, 
relevant results obtained in the RWA are plotted in Figs. [2](d)-(f ). Exact calculations in 
Fig-! a) show that for a held with / = 0 . 01 , magnitude of |a 0 (t )| 2 is decreased while that of 


11 







exact calculation 


RWA 



FIG. 2: (Color online) Time developments of |ao(t )| 2 (solid curves), |ai(f )| 2 (dashed curves) and 
1 03 (f) | 2 (bold solid curve) for (a) / = 0.01, (b) / = 0.02 and (c) / = 0.05 in exact calculations, and 
those for (d) / = 0.01, (e) / = 0.02 and (f) / = 0.05 in the RWA (g = 0.01 and u> = A 10 ), bottom 
curves in (a)-(f) expressing applied fields. 

|ai(t)| 2 is increased at t ~ 0 with small |« 3 (^) | 2 , which is similar to results in the RWA shown 
in Fig. 2(d). For / = 0.02, however, magnitude of |a 3 (f)| 2 in exact calculations becomes 
appreciable [Fig. [2j(b)] whereas it is vanishing in the RWA [Fig. |2](e)]. A comparison between 
Figs. [2j[c) and (f) show that the difference between an exact calculation and the RWA is 
evident for / = 0.05. 

Figures^ a), (b) and (c) express time dependences of (X 1 +X 2 ) for / = 0.01, 0.02 and 0.05, 
respectively, obtained by exact calculations (solid curves) and the RWA (dashed curves). We 
note that ()Ci + x 2 ) shows a complicated time dependence which arises from a superposition 
of multiple motions with frequencies of u and u ± as the RWA analysis shows. This 
analysis may be applied to case of / = 0.01 and / = 0.02. For a larger / = 0.05, however, 
(xi + X 2 ) of the RWA is very different from that of an exact calculation [Fig. 0c)]. 
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FIG. 3: (Color online) (a)-(c) {X 1 +X 2 ) for (a) / = 0.01, (b) / = 0.02 and (c) / = 0.05; (d)-(f) T(i) 
for (d) / = 0.01, (e) / = 0.02 and (f) / = 0.05; (g)-(i) C(t) for (g) / = 0.01, (h) / = 0.02 and (i) 
/ = 0.05, solid and dashed curves expressing results of exact and RWA calculations, respectively 
(g = 0.01 and u = A 10 ). Bottom curves in (a)-(i) denote applied fields. 

Figures |3](d), (e) and (f) show the correlation T(t) calculated for / = 0.01, 0.02 and 0.05, 
respectively, with g = 0.01. RWA analysis shows that T(f) 2 oscillates with frequencies of uj, 
and uj ± fh 

The concurrence C(t) calculated for / = 0.01, 0.02 and 0.05 are shown Figs. |3](g), (h) and 
(i), respectively. At t — 0, C'(O) = 0.1485. C(t) shows complicated time dependence because 
C(t) 2 includes superposed oscillations with frequencies of cu, 2u, Q, 20, uj ± 0, 2u ± fl and 
2(a; ± 0), as analyzed by the RWA which is expected to be valid for a small /. 

B. g dependence 

Next we report calculated results when the coupling g is varied. Figure |4](a), (b) and (c) 
show time developments of \a u \ 2 for g = 0.01, 0.1 and 0.2, respectively, obtained by exact 
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FIG. 4: (Color online) Time developments of |ao(f)| 2 (solid curves), |ai(f)| 2 (dashed curves) and 
|d 3 (t) | 2 (bold solid curve) for (a) g = 0.01, (b) g = 0.1, and (c) g = 0.2 in exact calculations, and 
those for (d) g = 0.01, (e) g = 0.1, and (f) g = 0.2 in the RWA (/ = 0.01 and u = Aio), bottom 
curves in (a)-(f) expressing applied fields. 

calculations with / = 0.01 and oj = A 10 (=0.07431, 0.02611 and 0.0140 for g = 0.01, 0.1 
and 0.2, respectively). For comparison, relevant results obtained in the RWA are plotted in 
Figs. |4](d)-(f ). A result of the RWA for g = 0.01 in Fig. |4](d) is in fairly good agreement 
with that of an exact calculation in Fig. |4](a). For g = 0.1, however, an agreement between 
exact and RWA calculations is not satisfactory [Fig. HR)]. For g = 0.2 result of the RWA 
is quite different from that of an exact calculation [Fig. 0C)]. 

The difference between results of | | 2 in an exact calculation and the RWA reflects on 

expectation values of (x\ + X 2 ) shown in Figs|5|a)-(f). We note that although expectation 
values in an exact calculation and the RWA are in fairly good agreement for g = 0.01, both 
the results are rather different for g — 0.1 and 0.2. 

Figs. [5j[d)- (f) and Figs. [5j[g)- (i) show the correlation and concurrence, respectively, 
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FIG. 5: (Color online) (a)-(c) {x\ + X 2 ) for (a) g = 0.01, (b) g = 0.1 and (c) g = 0.2; (d)-(f) T(i) 
for (d) g = 0.01, (e) g = 0.1, and (f) g = 0.2; (g)-(i) C(t) for (g) g = 0.01, (h) g = 0.1, and (i) 
g = 0.2, solid and dashed curves expressing results of exact and RWA calculations, respectively 
(/ = 0.01 and u = A10). Bottom curves in (a)-(i) denote applied fields. 


obtained by an exact calculation and the RWA for several g values with / = 0.01. We note 
in Figs. |5](g)-(i) that the concurrence becomes larger for larger g. In fact, the time-averaged 
(CRWA(t) 2 )av given by Eq. (83) is 0.383265, 0.634747 and 0.712556 for / = 0.01, 0.02 and 
0.05, respectively. 


C. c 0 dependence 

Although we have so far assumed that c 0 is equal to A 10 , its value is changed in this subsec¬ 
tion. Figures |6|a), (b) and (c) show exact calculations of time-dependent level populations 
\a u (t)\ 2 for cu/Aio = 0.8, 1.0 and 1.2, respectively, with f — g — 0.01. For comparison, 
relevant results in the RWA are plotted in Figs. [6](d)-(f ). In the resonant condition with 
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FIG. 6: (Color online) Time developments of |ao(t)| 2 (solid curves), |ai(t)| 2 (dashed curves) and 
|ci 3 (t) | 2 (bold solid curve) for (a) w/Aio = 0.8, (b) 1.0, and (c) 1.2 in exact calculations, and those 
for (d) w/Aio = 0.8, (e) 1.0, and (f) 1.2 in the RWA (/ = g = 0.01), bottom curves in (a)-(f) 
expressing applied fields. 


uj = A 10 , |a 0 (t)| 2 and |ai(t)| 2 oscillate between 0 and 1. It is not the case in the off-resonant 
condition with to ^ A 10 . We note in Figs. 10(d) and (f) that a result for cj/A 10 = 0.8 is 
the same as that for cu/Aio = 1.2 in the RWA. This is because 17 is expressed in term of 


(oj — Aio) 2 in Eq. (54). In exact calculations, however, a result for cj/Aio = 0.8 is different 
from that for cu/Aio = 1.2 as shown in Figs. |6](a) and (c). In particular, an exact calculation 
for cn/A 10 = 1.2 in Fig. |6](c) shows a peculiar time dependence, which is quite different from 
a result of the RWA in Fig. |6](f). 

Time dependences of {x\ + x?) for cu/Aio = 0.8, 1.0 and 1.2 are shown in Fig. (b) 
and (c) obtained by exact calculations (solid curves) and the RWA (dashed curves). 

Figures |7](d)-(f) and Figs. [T](g)-(i) show the correlation function and concurrence, respec- 
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FIG. 7: (Color online) (a)-(c) {x\ + X 2 ) for (a) u/ A 10 = 0.8, (b) 1.0, and (c) 1.2; (d)-(f) T(t) for 
(a) w/Aio = 0.8, (b) 1.0, and (c) 1.2; (g)-(i) C(t ) for (a) uj/ A 10 = 0.8, (b) 1.0, and (c) 1.2, solid 
and dashed curves expressing results of exact and RWA calculations, respectively (/ = 0.01 and 
g = 0.01). Bottom curves in (a)-(i) denote applied fields. 

tively, for co/A w = 0.8, 1.0 and 1.2 calculated by exact calculations (solid curves) and the 
RWA (dashed curves). Again results of T(t) and C{t) in the RWA for cu/Aio = 0.8 are the 
same as those for uj/A w = 1.2. Results of the RWA are not in good agreement with those of 
exact calculations except for the resonant condition of oj = A 10 with a small / and a weak 
9- 
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FIG. 8: (Color online) Magnitudes of wavefunctions \^{x\,X 2 ,t)\ 2 at (a) t = 0.0, (b) t = 100.0, 
(c) t = 200.0 and (d) t = 300.0 with / = 0.02, g = 0.01 and w = Aio for the initial wavepacket 


given by Eq. (85) 


IV. DISCUSSION 


A. Initial wavepacket state 


It would be interesting to study effects of fields applied to the initial wavepacket state 
given by 


°o(0) 


ai(0) 


1 


^2(0) — 03(0) — 0. 


(85) 


Figure [8]( a) shows the 3D plot of magnitude of |\P(xi, x 2 , £)| 2 at t — 0.0 as functions of 
X\ and X 2 with / = 0.02, g = 0.01 and c 0 = Aio, which initially has a main peak at 
(xi,X 2 ) = (1-23,1.23). Figures |8](b)-(d) will be explained shortly. 

Figures |9]( a), (b) and (c) show time developments of populations of \a u (t)\ 2 in levels v 
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exact calculation 
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FIG. 9: (Color online) Time developments of |ao(f)| 2 (solid curves), |ai(f)| 2 (dashed curves) and 
103 (t) | 2 (bold solid curve) for (a) / = 0.0, (b) / = 0.01 and (c) / = 0.02 in exact calculations, and 
those for (d) / = 0.0, (e) / = 0.01 and (f) / = 0.02 in the RWA; the initial wavepacket given by 


Eq. (85) is adopted with g = 0.01 and oj = Aiq. Bottom curves in (a)-(f) express applied fields. 


(= 0,1,3) for / = 0.0, 0.01 and 0.02, respectively with g = 0.01 and oj = A 10 (= 0.07431) 
obtained by exact calculations. Relevant results obtained in the RWA are plotted in Figs. 
|9](d)-(f ). For / = 0.0, the populations are time independent, |ao(t)| 2 = |ai(f)| 2 = 1/2 and 
|«2(^)| J = 103(i) | 2 = 0.0. For / = 0.01, |a 0 (t)| 2 is initially increased while |ao(t)| J is decreased 
from their initial values of 1/2 by an applied field. A result of the RWA for / = 0.01 in Fig. 
[9](e) is in fairly good agreement with that of an exact calculation in Fig. [9](b). However, Figs. 
[9](c) and (f) show that for / = 0.02 both results are different, in particular, a population of 
|<3-3 (^) | “ becomes appreciable in an exact calculation while it is vanishing in the RWA. 

3D plots of magnitudes of wavefunctions |d/(rci, X 2 , t)\ 2 at t — 100.0, 200.0 and 300.0 
are presented in Figs. |8](b), (c) and (d), respectively. At t = 100.0, the wavepacket has 
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FIG. 10: (Color online) (a)-(c) ( x\ + x 2 ) for (a) / = 0.0, (b) / = 0.01 and (c) / = 0.02; (d)-(f) 
F(t) for (d) / = 0.0, (e) / = 0.01 and (f) / = 0.02; (g)-(i) C(t) for (g) / = 0.0, (h) / = 0.01 and (i) 
/ = 0.02, solid and dashed curves expressing results of exact and RWA calculations, respectively, 


for the initial wavepacket given by Eq. (85) (g = 0.01 and u = Aio). Bottom curves in (a)-(i) 
denote applied fields. 


four small peaks at (xi,x 2 ) = (±1.23, ±1.23) and (xi,x 2 ) = (±1.23, =f 1.23) [Fig|8](b)]. At 
t = 200.0, it has again four small peaks at the same positions as at t = 100.0 [Fig|8|c)] 
although a peak at (aq, x 2 ) = (—1.23, —1.23) is the highest among the four. The wavepacket 
returns approximately to the initial state at t = 300.0 as shown in Fig. [8](d), which is similar 
to Fig. |8|a) at t = 0.0. 

Figures [lofa), (b) and (c) express time dependences of (xi+x 2 ) for / = 0.0, 0.01 and 0.02, 
respectively, obtained by exact calculations (solid curves) and the RWA (dashed curves). We 
note in Fig. [l0|(a) that ( X\ ± x 2 ) for / = 0.0 shows the simple sinusoidal time dependence, 
which express a tunneling of wavepacket between two minima of DW potential. When a 
field with / = 0.01 is applied, both results of (aq ± x 2 ) obtained by exact calculations and 
the RWA show complicated time dependences [Fig. [l0|(b)] . Although the both results are 
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in good agreement for / = 0.0 and 0.01, results of exact and the RWA calculations become 
different for a larger / = 0.02 as shown in Figs. [!o|c). We note that a tunneling becomes 
difficult by applied holds, which is nothing but a suppression of tunneling by coherent fields 

0 - 

Figures [To^d) , (e) and (f) show the time dependence of correlation T(t) for / = 0.0, 0.01 
and 0.02, respectively, with g = 0.01 obtained by an exact calculation (solid curves) and the 
RWA (dashed curves). r(t) for the adopted wavepacket with / = 0.0 is given by [12] 

r(f ) 2 = \[ l + cos(Ai 0 t)], (86) 


which leads to a sinusoidal oscillation with a period of 27 t/Ai 0 = 84.55. Figure |To|(e) shows 
that T(f) for / = 0.01 shows complicated time dependences and that a result of the RWA is 
in fairly good agreement with that of an exact calculation. For a larger / = 0.02, however, 
a result of the RWA becomes different from that of the exact calculation [Fig. [To|(f )]. 

Figures [lQ|(g) , (h) and (i) show the concurrence C(t ) for / = 0.0, 0.01 and 0.02, respec¬ 
tively, with g = 0.01 obtained by an exact calculation (solid curves) and the RWA (dashed 
curves). C(t) of the adopted wave packet for / = 0.0 is given by [12] 


C(t ) 2 = -[1 + sin 2 (26>) — 2 sin(20) cos(2A 10 f)], 


(87) 


which yields C(0) = 0.42577 at t = 0.0 for g = 0.01. C[t) for / = 0.0 oscillates with the 
period of vr/A 10 = 42.28 as shown in Fig. 10 g). Figures |!o|(h) and (i) show that for finite 
/ = 0.01 and 0.02, C{t ) has complicated time dependence, just as r(t). 


B. Step field 

Although sinusoidal periodic fields have been so far adopted, we may apply our method 
to any time-dependent field. As a typical example, we employ a step field given by 

F(t) = fQ(t ), (88) 

where 0(f) denotes the Heaviside function. 

We first try to obtain approximate analytical expressions of a„(t) (u — 0 — 3) for an 
applied step field with the use of the two-level approximation (TLA) where contributions 
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FIG. 11: (Color online) Time developments of |ao(t)| 2 (solid curves), |ai(f)| 2 (dashed curves) 
and |a- 3 (t)| 2 (bold solid curve) for the applied step field with (a) / = 0.01, (b) / = 0.02 and (c) 
/ = 0.05 in exact calculations, and those for (d) / = 0.01, (e) / = 0.02 and (f) / = 0.05 in the 
TLA (g = 0.01), bottom curves in (a)-(f) expressing applied step fields. 


only from two terms of a 0 (t) and Gq(f) are included. Equations (36)-(39) for t > 0 become 


daft 

in —— 
dt 

= —afe lAlot ai 1 

(89) 

dci\ 
in —— 
dt 

= —af e iAlot a 0 . 

(90) 

d(X 2 

da* 


in —— 
dt 

= ih—^ = o. 

dt 

(91) 
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FIG. 12: (Color online) (a)-(c) {x\ + X2) for (a) / = 0.01, (b) / = 0.02 and (c) / = 0.05; (d)-(f) 
r(f) for (d) / = 0.01, (e) / = 0.02 and (f) / = 0.05; (g)-(i) C(t) for (g) / = 0.01, (h) / = 0.02 and 
(i) / = 0.05, solid and dashed curves expressing results of exact and TLA calculations, respectively, 
(g = 0.01). Bottom curves in (a)-(i) denote applied step fields. 


Solutions of Eqs. (89)-(91) are given by 


a°(t) = [(A 10 -fI s )Ae - i ( Alo+n ^/ 2 + (A 10 + ^Se- i ( Alo -^) t / 2 ], 

a, (t) = A e ‘ (A “ -n ' ),/2 + B e i(Al " +n ' l</2 , 

a 2 (<) = 02(0), 

03 (t) = <13(0), 


(92) 

(93) 

(94) 

(95) 


with 

L! s = \J Af 0 + 4(af /h) 2 , (96) 

where A and B stand for integration constants. For assumed initial conditions ao(0) = 1 


23 










































































































(97) 


and ai(0) = 02(0) = <23(0) = 0 which yield A = —B = —af/hfl s , solutions are given by 
a 0 (t) = (^j [(A 10 + n 8 ) e -*(Aio-n.)‘/2 _ ( Al0 _ Q s ) e -i(A 1 o+n.)t/2] j 


ai(t) = - {w, 

a 2 (t ) = a 3 (t) = 0. 


3 *(Aio— n s )t/2_ i(Ai 0 +n ( 


)*)/ 2 ] ; 


(98) 

(99) 


Then |ao(t)| J and |ai(t)| J are expressed by superposed oscillations with frequencies of ± 


A 


lo- 


Figures nTl a), (b) and (c) show \a u \ 2 of exact calculations for step fields with / = 0.01, 


0.02 and 0.05, respectively, which are applied to the initial ground state given by Eq. (84); 
relevant results in the TLA are plotted in Figs. B d ). (e) and (f). When a step held is 
applied, \a u (t)\ 2 begin to oscillate. Equations (97) and (98) may explain time dependences of 
|ao(£)| J and |ao(£)| J for small / = 0.01 and 0.02. They are, however, not valid for / = 0.05, 
for which |a3(t)| 2 has an appreciable magnitude while it is assumed to be zero in the TLA. 


Time dependences of (aq +x 2 ), T(i) and C(t) are plotted in Figs. |l2[a)-(c), Figs. 12(d) 


(f) and Figs. 12(g)-(i), respectively, where solid and dashed curves denote results in exact 
calculations and in the TLA, respectively. Their time dependences for / = 0.01 and 0.02 


may be approximately elucidated in the TLA given by Eqs. (97) and (98), although they 
are not applicable to the case of / = 0.05. 


V. CONCLUSION 

We have studied effects of applied fields in quantum coupled DW system with the use 
of an exactly solvable Razavy’s potential [9]. From the Schroclinger equation for the driven 
DW system, we have obtained equations of motion for populations of the four levels. Model 
calculations of expectation values (xi+x 2 ), correlation T(f) and concurrence C(t) for applied 
sinusoidal fields show very complicated time dependence. Their time dependence may be 
analytically understood within the RWA in cases of a weak coupling and a small held in the 
near-resonance of 1 00 — A 10 | 00 . Otherwise, results of the RWA are not in good agreement 

with exact numerical calculations. It is indispensable to develop an analytical method going 
beyond the RWA for the DW system, just as for the TL model p7jRTT] . In the present study, 
we do not take into account environmental effects which are expected to play important 
roles in real DW systems. These are left as our future subjects. 
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Appendix: A. General driving fields 

We consider the input term which is more general than Eq. Q, as given by 

Hi = ~XiFi(t) — x 2 F 2 (t), (Al) 

F n {t) (n = 1,2) expressing a time-dependent field applied to the nth DW system. The 
energy matrix of the time-dependent total Hamiltonian H (= Hq + He + Hi) with basis 
states of <f>o, $1, $2 and $3 is given by 


H = 


Eo -f[F,(i) + F 2 (f)] -f[U(i) - F 2 (i)] 

-f[fi(«) + r 2 (i)] £1 0 

-finW-fiW] 0 e 2 

0 -f[*i(i) + F 2 (t)] f[nw - f 2 («>] 


0 \ 

-f[n(«) + r 2 (i)] 
f [nw - f 2 («)] 

) 


(A2) 


where a and [5 are given by Eqs. (29) and (30), respectively. Note that in the case of 


Fi(t) = F 2 (t) — F(t), Eq. (A2) becomes Eq. (28). 


With the use of Eqs. (34) and (A2), equations of motion for a^t) (p = 0 — 3) become 

(A3) 

” a 3, 


“IF = -|[U(i) + n(«)]e-‘ a '» , o 1 -2[F 1 ( t )-F 2 ( t )]e-<A»< 02 , 


= -|[n(() + F 2 («)]e ,a “‘a„-|[nM + r 2 («)]e- iAa “' 


<fi ‘E = -|[n(t) - F 2 (f)] e ia ”‘a„ + |[F,(i) - F 2 (i)] e- ia ”‘ a3 , 
th'E = -|[F 1 (() + F 2 (()]e ia »“a 1 + |[F 1 (()-F 2 («)]e‘ a31, a 2 . 


(A4) 

(AS) 

(A6) 


In the symmetric case of iq (£) = F 2 (f) = F(t), Eqs. (A3)-(A6) reduce to Eqs. (36)-(39). 
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FIG. 13: (Color online) Time developments of |ao(i)| 2 (solid curves), |ai(f)| 2 (dashed curves), 
|a 2 (t)| 2 (chain curves) and |a- 3 (i)| 2 (bold solid curve) with (a) / = 0.01, (b) / = 0.02 and (c) 
/ = 0.05 in exact calculations, bottom curves in (a)-(c) expressing applied fields of F\ (t) [F 2 (t) = 0]. 


In the anti-symmetric case of iq(f) = —F 2 (t) = F(t ), Eqs. (A3)-(A6) become 


ih —F = —f3F(t) e 
at 

ih ^ = o, 

dt 
da 


—i/±.2Qt 


« 2 ; 


ih —- = —/ 6F(t ) e* A2ot a 0 + aF(t)e 
dt 

ih = aF(t) e* As2t a 2 . 

LLL 


3 -?A 32 t 


a 3j 


(AT) 

(AS) 

(A9) 

(A10) 


for which we obtain a\(t) = ai(0). 

When helds are applied only to the first DW system with iq (t) ^ 0 and F 2 (t) = 0, Eqs. 
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FIG. 14: (Color online) (a)-(c): Time dependences of (xi) (solid curves) and ( 2 : 2 ) (dashed curves) 
for (a) / = 0.01, (b) 0.02 and (c) 0.05; (d)-(f): T(t) for (d) / = 0.01, (e) 0.02 and (f) 0.05; (g)-(i): 
C(t) for (g) / = 0.01, (h) 0.02 and (i) 0.05 (g = 0.01). Bottom curves in (a)-(i) denote applied 
fields of F\ (t) [F 2 (t) = 0]. 


(A3)-(A6) yield 


dao 
in —— 
dt 

dai 
in —— 
dt 

da,2 
in —— 
dt 

da 3 
in —— 
dt 


For applied fields of F { (t) = 
and a 3 (t) are negligible as will 



(All) 

aFl ^ e iAlot a 0 l5Flif) e- iA ^a 3 , 

2 2 

(A12) 

Pnm e .A 2 „i 0o + e -iA,al 03j 

(A13) 

w(i) + aFl(t) 

2 2 

(A14) 

/ sin cut and F 2 (f) = 0 with small / and weak g where a 2 (t) 
be shown shortly, we may take into account only a 0 (t) and 
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ai(t). Then Eqs. (All)-(A14) within the RWA become 


dao 
dt 
dai 
dt 


af_ 

4 h 


e i(u)-A 10 )t 


af_ 
4 h 


e -i(u-A 10 )t a()) 


da 2 del 3 

dt dt 


(A15) 

(A16) 

(A17) 


which are equivalent to Eqs. (45 )-(4T) for the symmetric driving fields if we read //2 —> /. 
Figures |l3[a)-(c) show time developments of \a u (t)\ 1 2 3 4 5 6 (v — 0 — 3) calculated for / = 0.01, 


0.02 and 0.05, respectively, with g = 0.01 and u = Ai 0 when fields of F\ (t) = f smut and 
F 2 (t) = 0 are applied to the system with an initial stable state given by Eq. ( [84] ) . When a 
field with / = 0.01 is applied, |a 0 (^)| 2 is decreased from unity while |ai(£)| J begins to increase 
at f ~ 0.0 with vanishingly small |cz. 2 (t) | J and |« 3 (^) | 2 • For larger fields with / = 0.02 and 


0.05, however, magnitudes of \a 2 {t)\ 2 and |a 3 (t )| 2 become appreciable. Results in Fig. 13)b) 
with / = 0.02 are not dissimilar to those in Fig. |2](a) for the symmetrical field with / = 0.01, 
as mentioned above. 

Time dependences of (aq) (solid curves) and (x 2 ) (dashed curves) for / = 0.01, 0.02 and 
0.05 are plotted in Figs. [I4|(a), (b) and (c), respectively, where (aq) is different from (x 2 )\ 
the latter seems a little delayed than the former. Note that (x 2 ) = 0.0 for g = 0.0 for which 
the two DW systems are decoupled (related results not shown). 


Relevant results of T(f) for / = 0.01, 0.02 and 0.05 are shown in Figs. 14'd), (e) and (f), 


respectively, while C(t) are plotted in Figs. 14 'g), (h) and (i). They show very complicated 
time dependence. It is impossible to theoretically elucidate time dependences of \a u (t)\ 2 , 
(x n ) (n = 1 , 2 ), T(f) and C(t) except for cases with very small / and weak g. 
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